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Abstract 

We study magnetic properties of itinerant quantum magnetic particles described by a generalized 
Hubbard model with large spin (S > 1/2) which may be realized in optical lattices of laser-cooled 
atom systems. In fermion systems (half-integer spins), an extended form of Nagaoka ferromag- 
netism may be realized. However, as novel aspects of the large spin cases, we found that the 
condition on the lattice connectivity is more stringent than in the case of 5 = 1/2 particles, and 
that the system shows a peculiar degenerate structure of the ground state in which the ferromag- 
netic state is included. In contrast, it turns out that the ground state of itinerant bosonic systems 
(integer spins) has a degenerate structure similar to that of fermion system with S* > 1/2 regardless 
of the shape, connectivity or filling of the lattice, and that the state with the maximum total spin 
is always one of the ground states. 

PACS numbers: 75.10.Jm,75.45.+j,75.75.+a,75.30.Gw,75.40.Mg,75.50.Xx 
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I. INTRODUCTION 



The origin of magnetism is attributed to the quantum mechanical interaction of particle 
swhich carry spin. In so-called localized spin systems, the Pauli principle plays an essential 
role and the magnetic interaction is expressed as the Heisenberg interaction^, the exchange 
integral between atoms being the dominant term. Not only the two-spin interaction but also 
multi-spin interactions may contribute to give rise to exotic magnetic states. In particular, 
various effects of the multi-spin interaction have been reported in ^He^. For itinerant electron 
systems, the origin of the magnetic order, and of ferromagnetic order in particular, has 
been studied extensively. Itinerant electron systems are often described by tight binding 
models such as the Hubbard model^. Nagaoka pointed out that the ground state of the 
Hubbard model may be a ferromagnetic state (Nagaoka ferromagnetic state) if the number 
of electrons is reduced by one from half-filling^i^, where half-filling means that the number 
of particles is the same as the number of the lattice sites. The Nagaoka ferromagnetic state 
is established if the so-called "connectivity condition" on the lattice is satisfied^i^. In the 
Nagaoka ferromagnetic state, the energy of the system is minimized when the total spin 
of the system takes the maximum value. On the other hand, the ground state of the half 
filled system is a singlet state, that is its total spin is zero. Elsewhere, we have studied how 
the magnetic state changes between these two states when an electron is removed from the 
system and demonstrated an adiabatic change between these states^. 

Magnetism has primarily been studied in electron systems in which the spin S = 1/2. 
However, recently developments in the field of laser-cooled atomic systems have opened 
new possibilities to realize artificial tight-binding quantum system such as the Hubbard 
model^iii. In contrast to electron systems, in optical lattices the spin of trapped atoms is 
not necessarily S* = 1/2 but can take larger half- integer or integer values. In the latter case, 
the system contains bosons, not fermions, and it is therefore of interest to study itinerant 
magnetism of systems with S* > 1/2. The present paper presents the results of such systems. 
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II. MODEL 



We consider a tight binding model of the form 

N 

<ij>,M i=l 

where c] and are the annihilation and creation operators of a particle (fermion or 
boson depending on its spin) of the magnetization M at a site i, respectively, rii^M is the 
number operator 

for the particle of spin S, M = —S,—S + l,...,^, and U^rii^pf) represents the on-site 
repulsive interaction. For the system with S = 1/2, we take for Uijii^M) the standard form 

U{ni^M) = Uoni_i/2ni^i/2, (3) 

and for larger 5* > 1 /2 we extend it as 

t/(n,,M) = ^UoNiiNi - 1), (4) 

where Ni = '^j^j rii^M and Uq is assumed to be positive. Obviously, the interaction Uijii^M) 
increases the energy whenever a site is occupied by more than one particle. 

In order to study magnetic properties of the system, we introduce spin operators for the 
magnetization S = (S'f , Sf , S-) of the particles, where 

s 

M=-S 

Hereafter, we use the operators = + iSf and S~ = — iSj which are expressed in 
terms of c\^i^j and c^^j^: 

St = Y ViS-M){S + M + l)cl^^,c,^^ (6) 

M 

and 

S- = J2 V(5 + M)(^-M + l)4^,_,c,^, (7) 

M 

where 5* is the total spin of each particle, and M is the ^-component of the magnetization. 
To discuss the magnetic properties of the states, we adopt the usual notation 15*, M) where 



S is the total spin 5* and M is the magnetization. The action of the operators Eqs. ([6]) and 
([7]) on the state IS", M) is given by the relations 



S^\S,M) = ^/iS-M)iS + M+l)\S,M + l) (8) 

and 

S-\S,M) = ^{S + M){S - M +1)\S,M -1). (9) 

In order to explicitly compute matrix elements of the operators, it is convenient to intro- 
duce orthornormal basis states. In the case of fermion systems, we adopt the form 

i^) = 4m---4m'|o)' (10) 

where i > j, and M > M' ii i = j. With this definition, the operations Eqs. (El) and (I7j) do 
not change the order of creation operators in Eq. ffTOl) . In the case of bosons, more than two 
particles with the same M can occupy the same site and the normalized basis states take 
the form 

I*) = / r = \ni,M)---\njM'). (11) 

For boson systems, the order of i and j and M and M' is not relevant. 
The total spin of the whole system is given by 

52 = (5^)2 + {Syf + (5^)2 = ^ + (12) 

where = 5f , 5^ = ^. Sf, and 5^ = ^. 5f , and 5± = 5f ± z^f . We denote 
the value of the total spin of the system by Stot, i-e., 5'tot('S'tot + 1) = {S'^). 

III. CONSERVATION OF NUMBER OF PARTICLES OF DIFFERENT SPINS 
AND GROUND-STATE DEGENERACY 

It should be noted that the Hubbard Hamiltonian conserves number of particles 

N 

nM = ^ni^M (13) 

i=l 

for each M. To describe the set of hm of states, it is convenient to introduce the notation 

{um} = {ns,ns-i,- ■ ■ ,n-s), (14) 



where Y1m=-s^^ ~ important to note that except for S = 1/2, the operator 

S^Sj' + S^S^ changes the set {um}- For example, if S* = 1 apphcation of S^Sj' + 'S'j~5'^ to 
states in the set (ni = 0, rio = 2, n„i = 0), created states in the set (ni = 1, no = 0, n„i = 1). 
Thus, except for S = 1/2, the matrix element {{n\f}\S^Sj' + Sj^\{nM}) can be nonzero 
even if {um} 7^ {?^a/}- 

Thus, even though the Hamiltonian Ti. has the SU(2) symmetry, conserves the total spin, 
and therefore commutes with S*^, the operation of S^Sj' + S^S^^ on each state with given 
{um} must be treated carefully. 

Let us denote by \G, M) the ground state in the space with the magnetization M . Because 
the Hamiltonian Ti and S"^ commute with each other, S^\G^M) is also a ground state 
because 

'HS~\G, M) = S~n\G, M) = EgS~\G, M). (15) 
However, the set of numbers {um} is not necessarily conserved, that is 

S-n\G,M, {nM])=EGS-\G,M, {um}) 

= EcJ2 «Km}I^' ^ - 1' {</))' (16) 

where \G, M, {um}) denotes one of basis states with fixed {um} that contribute to the 
expansion of \G, M) in terms of basis states. The energy of this state is given by 

{G,M,{nM}\S+nS-\G, M, {um}) 



G 



{G,M,{nM}\S+S-\G,M,{nM}) 



Because 



^ Ekj Hn'J'jG, M - 1, {n\,}\n\G, M - 1, {n'^}) 

EkjI«k,}P(g,m-i,{<,}|g,m-i,k,}) ■ ^ ' 

{G,M-l,{n'j,,}\n\G,M-l,{n'M}) ^ p 
(G,M-1,K,}|G,M-1,K,}) - ^ ^ 

in order to satisfy the relation f|T7j) . the state in each set of numbers {n\j} must be the 

ground state, i.e. 

{G,M-l,{n\,}\n\G,M-l,{n',,}) 
(G,M-l,K/}|G,M-l,{n^}) ^' ^ ^ 

and thus 

n\G, M-l, {n',j}) = Eg\G, M - 1, {n\j}). (20) 
This shows that the same ground state energy is found for all sets {n'^j} except if a^n'j^^} = 0. 



5 



(a) (b) 
FIG. 1: Lattices: (a) five sites and six bonds; (b) lattice (a) with two additional bonds. 
IV. GROUND STATE OF A FERMION SYSTEM WITH 5 = | 
A. Subspaces due to conservation of particle number of each M 

Let us now study the dependence of the ground state energy on Uq- We consider a system 
with four particles with S = 3/2 on the 5-site lattice depicted in Fig.HJ^a). Figure [2](a) shows 
the ground state energies for the sets {?t.a/}=('^3/2 5 ^^1/25 ^^-1/25 ''^-3/2) = (4, 0, 0, 0), (3, 1, 0, 0), 
(2, 2, 0, 0), (2, 1, 1, 0), (1, 1, 1, 1), corresponding to systems in which the number of different 
states is 5, 50, 100, 250 and 625, respectively. The Hamiltonian does not depend on the 
values of M , and thus the energy level structure is same for the cases with the same set of 
numbers, e.g. for (3,1,0,0),(3,0,1,0),- ■ ■ , (0,0,1,3). For small value of Uq, the ground state 
energies are all different. As Uq increases, the ground state energies of (3,1,0,0), (2,2,0,0) 
and (2,1,1,0) become degenerate with that of (4,0,0,0). This fact indicates that the lowest 
energy state in these sets has the total spin ^max- However, we find that the lowest energy 
state in the set (1,1,1,1) is always lower than that of (4,0,0,0). This means that even at large 
Uq the ground state has total spin S < S^ax- Therefore, the Nagaoka-ferromagnetic state, 
i.e., the state of the maximum total spin, is not realized as the ground state in the present 
case even at large values of Uq. 

In general, when there are multiple conserved quantities in a system that do not commute 
with each other, each energy state of the system is usually degenerate as shown in Eq. (120P ^. 
In the present model, the total magnetization and the set {um} are conserved. However, the 
ground state in Fig. 2(a) is not degenerate, which gives an exception to the above general 
property. 

Figure 3 shows the Uq dependence of the total spin of the lowest-energy state for the 



6 



sets {um} on the lattice of Fig. [2](a). The total spin of states in the subspace (4,0,0,0) 
is Sma,x = (3/2) X 4 = 6. For small values of U, the total spin of states in the subspace 
(3,1,0,0) is less than S'max, but it becomes S'max when the ground state energy becomes 
degenerate to that of the (4,0,0,0) subspace. The ground state energies of (2,2,0,0) and 
(2,1,1,0) become degenerate to that of (4,0,0,0) at certain values of Uq. However, the total 
spin of the lowest-energy state of these sets does not reach S'max- This fact agrees with the 
earlier argument that if S~\G,M) consists of more than one set of {riA/j's, neither of these 
sets yields an eigenstate of the total spin although each of them have the same ground state 
energy. Therefore, the expectation value of the total spin is not necessarily an integer. 

B. Connectivity condition for the system with S = | 

Note that the lattice shown in Fig. [D^a) satisfies the so-called connectivity condition for 
the case of 5 = 1/2 and the Nagaoka-ferromagnetic state is realized at large Uq. However, 
for S = 3/2, this state is not realized. The reason is that the connectivity condition for 
four S = 3/2 particles with different magnetization, i.e., M = —3/2,-1/2,1/2,3/2 is not 
satisfied in the lattice of Fig. [It^a). Indeed, the configuration depicted in Fig. |H^a) can not 
be changed into that of Fig. 111(b). 

In order to re-establish the connectivity condition we may add bonds and construct 
the lattice depicted in Fig. [H^b). This lattice has loops with an odd number of bonds and 
therefore we need to use a negative value of t. Then, we find that the Nagaoka-ferromagnetic 
state is realized, as is clear from Fig. [2t^b). Summarizing, we have shown that S = 3/2 
itinerant particle systems may exhibit Nagaoka-ferromagnetism although it is more difficult 
to satisfy the connectivity condition. 

C. Structure of the ground state 

We found that the ground state in each subspace for a set {um} is unique but in each 
subspace we have one state with the same ground-state energy. In Table HI we list all the sets 
{um}- The structure indicates that there is a state with S'tot = 6,4, 2 and 0. In the case of 
S = 1/2, there is a one-to-one correspondence between the magnetization M and {n+,n_}, 
as is clear from Table [B Thus, when we create states by applying S~ from the all-up state. 
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FIG. 2: Energies of the ground state for several sets of (^3/21 ?^i/2i ''^-1/21 "'^-3/2)- (a) A system with 
four particles on the lattice shown in Fig. H^a) . The solid line denotes the ground state energy for 
(4,0,0,0), circles: (3,1,0,0), squares: (2,2,0,0), triangles: (2,1,1,0), and bullets: (1,1,1,1). (b) Same 
as (a) but for the lattice shown in Fig. [IJb). 

we have only the state with Stot = 2 for the four spin system. In contrast, in the case of 
S = 3/2, we create different sets by applying S~. As we mentioned in the previous section, 
the created states are degenerate as ground states. Therefore, the structure displayed in 
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FIG. 3: The total spin S'tot of the lowest-energy state for the cases shown in FiglJl^a). The legend 
is the same as in Fig. 2. 




(a) 



(b) 



FIG. 4: Nearest-neighbor hopping cannot change configuration (a) into (b) without creating double 
occupancy. 

Table [H is intrinsic for systems with S = 3/2, 5/2, . . .. 

V. GROUND STATE OF A BOSON SYSTEM WITH S = 1 

Next, we study a system of S" = 1 particles, namely a boson system. Figure [5] shows the 
ground state energies for the sets {um} = {ni,no,n^i)= (4,0,0), (0,4,0), (3,1,0), (3,0,1), 
(2, 2, 0), (2, 1,1), (1, 2, 1), corresponding to systems in which the number of different states 
is 70, 70, 175, 175, 225, 375, and 375, respectively. Surprisingly, we find that the energies are 
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TABLE I: Sets of particles of different M{> 0) {um} for the cases oi S = 3/2 and 1/2. For 
S = 3/2, we list the sets only for positive M. 



5 = 3/2 




M 


6 5 4 3 2 1 


W/} 


(4,0,0,0) (3,1,0,0) (3,0,1,0) (3,0,0,1) (0,4,0,0) (0,3,1,0) (1,0,3,0) 
(2,2,0,0) (1,3,0,0) (2,1,0,1) (2,0,1,1) (0,3,0,1) 
(2,1,1,0) (1,2,1,0) (1,2,0,1) (2,0,0,2) 

\ 3 1 1 / \ 1 1 1 / \ 7 7 7 / \ 7 7 7 / 

(2,0,2,0) (1,1,2,0) (0,2,2,0) 

(1,1,1,1) 


S= 1/2 




M 


2 10-12 


{nu} 


(4,0) (3,1) (2,2) (1,3) (0,4) 



TABLE IL Sets of particles of different M(> 0) {nu} for the case of S" = 1. 



S = 1 




M 


4 3 2 1 


{nAi} 


(4,0,0) (3,1,0) (3,0,1) (2,1,1) (0,4,0) 




(2,2,0) (1,3,0) (2,0,2) 




(1,2,1) 



the same for all the cases regardless of Uq. The total spin of the different ground states are 
also plotted in Fig. [51 Thus, the magnetic properties of the boson system are very different 
from that of the fermion system. The total spins are less than the maximum value (S'tot=4) 
except for (4,0,0) and (3,1,0) which reflects the fact that the eigenstate of maximum total 
spin contains more than two sets of {um}, as in the case of the S = 3/2 system. 

We find a ground state in each set as in the case of S* = 3/2. We list the sets in Table ITTl 
Here we find again a degenerate ground states with ^tot = 4, 2 and 0, which is an intrinsic 
property of systems with S = 1,2,.... 

Moreover, we find that even in the half-filled case, in the boson system the total spin 
takes the maximum value. For instance, for a system of four atoms on a simple square 
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lattice (corresponding to the half-filled case), the ground state is the same for all sets 
(results not shown). 

This property of boson systems can be understood as follows. Consider the subspace for 
a fixed set {um}- All the off-diagonal matrix elements of the Hamiltonian 7i are —t or 0. By 
substracting an appropriate multiple of the unit matrix, also the diagonal elements can be 
made negative or zero. Thus, all the elements of the shifted matrix H are either negative or 
zero. In our model, all the sites are connected by bonds and therefore there exist a number 
n > such that all the elements of H^"' are positive. Then, the Perron-Frobenius theorem 
tells us that there exist a unique eigenstate of H^^ with an eigenvalue that is larger than 
the absolute value of all other eigenvalues. This unique eigenstate is therefore the ground 
state of H, is totally symmetric with positive coefficients, and contains the state of the 
maximum total spin. However, as we mentioned above, this ground state has an intrinsic 
degeneracy with respect to subspaces that have different {um} because we always start from 
the state with all spins maximum and let S~ create different sets of {n^j}. The states that 
are generated in this manner have all positive coefficients and are therefore ground states 
too. Cleary, this property does not depend on the connectivity of the lattice or the value of 
Uo. 

In contrast, for fermion systems, it is in general impossible to transform H such that all 
elements have the same sign but in those cases for which such a transformation exist, which 
is precisely the condition of Nagaoka ferromagnetism, we can apply the same arguments as 
in the boson case to prove that the ground state of H is totally symmetric with positive 
coefficients and contains the state of the maximum total spin. 

VI. SUMMARY AND DISCUSSION 

We have studied the magnetic properties of the ground state of itinerant systems with 
S > 1. We found that fermion systems (5* = 3/2) support an extended form of Nagaoka 
ferromagnetism but that the connectivity condition is more difficult to satisfy because of 
the presence of particles with different magnetization. When the maximum ^tot state is the 
ground state, the system has a degenerate ground state in each set {um} listed in Table [H 
Thus the ground state manifold consists of not only the state of the maximum Stot but also 
of states with smaller values of Stot- This degenerate structure is intrinsic for the systems 
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FIG. 5: Ground state energies E and the corresponding total spill 5'tot ^ 

function of Uq for four 

S* = 1 bosons and various (ni,no,ri-i) on the lattice shown in Fig. [T|^a). The solid line denotes 
data for (4,0,0), solid triangles: (3,1,0), open triangles: (3,0,1), open squares: (2,2,0), reversed 
triangles: (2,1,1), crosses: (1,2,1), and bullets: (0,4,0). 

with S > 1/2 where several sets of {um} exist for a given value of M. 

For boson systems (S* = 1) we found that there is a ground state in each set of {um} 
and that the same type of degenerate ground-state structure appears as the one found for 
5 = 3/2, but in this case regardless of the value of Uq and of the shape of the lattice, which is 
an intrinsic property of the bosonic case. Thus, we conclude that boson itinerant magnetic 
systems always have a state with the maximum total spin belonging to the manifold of 
ground states. This property follows quite naturally from the fact that the Hamiltonian of 
the boson system can be transformed such that all elements have the same sign, implying 
that the ground state is fully symmetrized. In contrast, for fermion systems we need an 
additional condition, the condition for the Nagaoka ferromagetism, for the ground state to 
have maximum total spin. 

We expect that these properties will be confirmed in real experimental systems. 
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